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Abstract. The phase transition to a Bose-Einstein condensate is unusual in that it 

' ' is not necessarily driven by inter-particle interactions but can occur in an ideal gas 

rrt as a result of a purely statistical saturation of excited states. However, interactions 

(JJQ are necessary for any system to reach thermal equilibrium and so are required for 

•j condensation to occur in finite time. In this Chapter we review the role of inter- 

r^ actions in Bose-Einstein condensation, covering both theory and experiment. We 

C^ focus on measurements performed on harmonically trapped ultracold atomic gases, 

3 but also discuss how these results relate to the uniform-system case, which is more 

WJ ' theoretically studied and also more relevant for other experimental systems. 

We first consider interaction strengths for which the system can be considered 
sufficiently close to equilibrium to measure thermodynamic behaviour. In particu- 

Ch lar we discuss the effects of interactions both on the mechanism of condensation 

. (namely the saturation of the excited states) and on the critical temperature at 

^ which condensation occurs. We then discuss in more detail the conditions for the 

^ equilibrium thermodynamic measurements to be possible, and the non-equilibrium 

^ phenomena that occur when these conditions are controUably violated by tuning the 

' ' strength of interactions in the gas. 
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1 Introduction 

Virtually all thermodynamic phase transitions are driven by interactions be- 
tween particles, which promote symmetry breaking into an ordered state. The 
phase transition comes about as a result of the competition between the en- 
ergy, which favours the ordered state, and the entropy, which favours the dis- 
ordered state. In contrast, Bose-Einstein condensation (BEC) is a purely sta- 
tistical phase transition, which at least in principle should not rely on interac- 
tions. The transition is instead a direct consequence of the finite-temperature 
saturation of the number of particles in the excited states of the system [1-4] . 
While this statistical argument does not explicitly invoke interactions between 
the particles, it does assume that the gas is in thermal equilibrium, which is 
impossible to attain in a completely noninteracting system "'^. This makes it 
challenging to experimentally observe ideal-gas behaviour and disentangle the 
role of interactions on the thermodynamics and dynamics of condensation. 

In this Chapter we review our recent experiments on this topic [6-8] , per- 
formed with an ultracold Bose gas of ^^K atoms with tuneable interactions. 
We were able to identify the interaction regime in which the gas may be con- 
sidered to be in thermal equilibrium and also to extrapolate our results to the 
noninteracting limit where direct equilibrium measurements are not possible. 
This allowed us to verify the statistical-saturation BEC mechanism in the 
noninteracting limit, and to accurately determine the deviations from ideal- 
gas behaviour due to interactions; these are seen both in the non-saturation 
of the excited states and in the shift of the critical point. Before presenting 
the experimental results we briefly review some background theory that will 
be useful for our discussion. 

1.1 Noninteracting Bosons 

We start by considering an ideal, noninteracting Bose gas. We first derive the 
key results for a uniform system, which we then apply to the trapped gas 
using the local density approximation (LDA). This "local" approach will be 
useful later when wc consider the effects of interactions, and in particular for 
the comparison of a uniform Bose gas with one that is harmonically trapped. 
The equilibrium momentum distribution of noninteracting bosons with 
mass m at a temperature T is given by the Bose distribution function 

where p is the momentum and /i < the chemical potential. The total particle 
density n can be found by integrating over all momentum states: 



^ In the recently observed Bose-Einstein condensation of a photon gas [5], there 
is no direct interaction between the light particles. However the interaction with 
the material environment, which ensures thermalisation, leads to a second-order 
interaction between the photons. 
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where 33/2(2;) — X]fe°=i ^'^/^'^ i^ ^ polylogarithni function and A = 
[2-Kh'^ / (mkBT)]^^^ is the thermal wavelength. We can re-express this result 
in terms of the phase space density D as 

D^nX^^gy.ic'^/'^-^). (3) 

Eq. (3) shows that there is a maximum value that D can take. This critical 
value is reached when /i = and is given by Dc = 33/2(1) — C(3/2) « 2.612 
(where C is the Riemann function). At a given temperature this corresponds 
to a maximum density. If this density is reached all the excited states saturate 
and any additional particles must accumulate in the ground state, forming a 
Bose-Einstein condensate^. At a given density n the BEC transition temper- 
ature is given by 

where the superscript ° refers to the fact this is an ideal gas result. 

For a gas in a potential V{r) we may apply LDA to Eq. (2). This amounts 
to having a local chemical potential 

fiir) = ^ - V{r) . (5) 

Specifically for a harmonic trap, V{r) = '^{l/2)mu!frf, where cui (with i — 
1,2,3) are the trapping frequencies along three spatial dimensions. 
The local density is then 

,, g3/2(e(^-^(-))/^-^) ,^, 

n{r) = -^ Y^ , (6) 

and the local phase space density D{r) = n{r)X^. The total number of particles 
in the excited states can be found by integrating over all space: 

, f 53/2(e(^-^("-'/'^-^) 
N' = ^-^^ '- dr . (7) 

The critical point for a trapped gas is the point at which the maximal 
local D reaches the critical value of ^(3/2). For a fixed T it makes sense to 
define the critical point in terms of the critical total particle number Nc- For 
a harmonic trap, with the geometric mean of the three trapping frequencies 
Q, the integral in Eq. (7) for /i = gives 



^ The singular ground-state contribution to the total density is implicitly excluded 
from the integral in Eq. (2). As /i approaches zero from below the ground state 
occupation can become arbitrarily large, as can be seen by inspecting Eq. (1). 
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Fig. 1. Ideal Bose gas condensation. Number of thermal atoms A'^' (black line) and 
number of condensed atoms A^'o (grey line) are plotted versus the total atom number 
A'^, at a fixed temperature. As atoms are added to trap N' = N and A^o ~ until the 
critical atom number A'^ is reached. At this point the excited states of the system 
saturate, and for N > N^ we have A^' = N° and No = N - N^^ . 
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where C(3) « 1.202. The equivalent expression for the transition temperature 
at a fixed particle number is given by^ 
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The ideal-gas picture of Bose-Einstein condensation driven by the purely 
statistical saturation of the excited states is simply summarised graphically 
in Fig. 1. Here we plot the number of atoms in the excited states, N' , and 
in the condensate, A^O) as the total atom number N is increased at constant 
temperature. For N < N^ no condensate is present and N' = N. However for 
N > iV° the thermal component is saturated at N' — N^ and the number 
of condensed atoms is simply given by A^o ~ N — N^. In sections 3 and 4 we 
will examine the effects of interactions both on the saturation of the excited 
states and on the value of the critical atom number. 



1.2 Interacting Bosons 

The dominant effects of interactions on Bose-Einstein condensation are quite 
different in a uniform system and in the experimentally pertinent case of a 



■^ Finite-size corrections slightly reduce the ideal-gas critical temperature, by 
ksATJ? = — ("(2)/(2(^(3)) fiojm ~ — 0.684/i(x)m, whore cjrn is the algebraic mean 



of the trapping frequencies [9] . 
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harmonically trapped atomic gas. This complex problem has a long history 
and for reviews we refer the reader to, for example, [3,4,9, 10]. Here we just 
briefly introduce some key points. 

The simplest theoretical framework to address the effects of interactions 
in a Bose gas is the Hartree-Fock approximation [9]. In this mean- field (MF) 
model one treats the thermal atoms as a "noninteracting" gas of density n'(r) 
that experiences a self-consistent MF interaction potential g[2no(r) -I- 2n'(r)], 
where g = Airh^a/m, a is the s-wave scattering length, and rio(r) the conden- 
sate density. We can then define an effective total potential 

Kfl-(r) - V{v) + 2.g[no(r) + n'(r)] , (10) 

and apply the LDA by replacing V{y) with V"off(r) in Eq. (6). Meanwhile the 
condensed atoms feel an interaction potential g[nQ{r) -I- 2n'(r)], where the 
factor of two difference in the condensate self-interaction comes about due to 
the lack of the exchange interaction term for particles in the same state"*. 

In a uniform system, the MF potential gives just a spatially uniform en- 
ergy offset and the most interesting effects arise due to beyond-MF quantum 
correlations. 

On the other hand, in a harmonically trapped gas (with repulsive interac- 
tions) the inhomogeneous density results in a mean-field repulsion of atoms 
from the central high-density region. This geometrical effect often dominates 
and makes it harder to experimentally observe the more interesting beyond- 
MF physics. 

1.3 Chapter outline 

In section 2 we bricfiy outline our experimental procedure for performing preci- 
sion measurements of the effects of interactions on Bose-Einstcin condensation 
of an atomic gas. 

In sections 3 and 4 we discuss the effects of interactions on the thermody- 
namics of a Bose gas with tuneable interactions. Here the range of interaction 
strengths we explore experimentally is such that the gas can always be as- 
sumed to be in thermal equilibrium. In section 3 we scrutinise the concept 
of saturation as the driving mechanism for Bose-Einstein condensation, and 
in section 4 we focus on the interaction shift of the critical point for con- 
densation. We compare the experimental results to both MF and beyond-MF 
theories and discuss how they relate to the case of a uniform Bose gas. 

In section 5 we discuss the conditions required for equilibrium measure- 
ments, and non-equilibrium effects that are observed when they are violated. 



■* This approach does not take into account the modification of the excitation spec- 
trum due to the presence of the condensate, which is included in more elaborate 
MF theories such as those of Bogoliubov [11] and Popov [12] (see also [9]). How- 
ever, it is often sufficient to give the correct leading order MF results. 
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2 Precision measurements on a Bose gas with tuneable 
interactions 

All the results presented here were obtained by performing conceptually simple 
experiments which are close in spirit to the ideal-gas theoretical plots of Fig. 1, 
but are performed for various strengths of repulsive interatomic interactions, 
characterised by the positive s-wave scattering length a. 

Our experiments start with a partially condensed gas of '^^K atoms, pro- 
duced in an optical dipolc trap [13] with u}/2tt varying between 60 and 80 Hz 
for data taken at different temperatures. The strength of interactions in the 
gas can be tuned by applying a uniform external magnetic field in the vicinity 
of a Feshbach scattering resonance centred at 402.5 G [14]. 

In each experimental series we fix a by choosing the value of the Feshbach 
field, and keep the temperature constant by fixing the depth of the optical 
trap. The total number of trapped atoms A^ is then varied by holding the gas 
in the trap for a variable time t, up to tens of seconds. During this time N, 
initially larger than the critical value Nc, slowly decays and eventually drops 
below Nc- Meanwhile elastic collisions between the atoms act to redistribute 
the particles between the condensed and thermal components of the gas. 

In each experimental run within a given series, corresponding to a par- 
ticular hold-time t, the thermal atom number N' and the condensate atom 
number Nq are extracted from fits to the absorption images of the gas after 
18 — 20ms of free time-of-fiight (TOF) expansion from the trap [15, 16]. The 
interactions are rapidly turned off at the beginning of TOF, by tuning the 
Feshbach field to the a — point. This minimises the condensate expansion 
and allows us to home in on the critical point by reliably measuring condensed 
fractions as low as 0.1 %. 

For accurate measurements of the small interaction shift of the critical 
point it is particularly important to minimise various systematic errors, e.g. 
due to finite-size effects [9], uncertainties in the absolute calibration of N 
and w, and small anharmonic corrections to the trapping potential [13]. We 
achieve this by performing "differential measurements", always concurrently 
running two experimental series which are identical in every respect except 
for the choice of the scattering length [7] . 



3 Non-saturation of the excited states 

In this section, we focus on the concept of the saturation of the excited states 
as the underlying mechanism driving the BEG transition. In superfiuid *He, 
which is conceptually associated with BEG, strong interactions preclude direct 
observation of purely statistical effects postulated by Einstein for an ideal gas. 
On the other hand it is generally accepted that a close-to-tcxtbook BEG is 
observed in the weakly interacting atomic gases. Therefore, one might expect 
that the saturation inequality N' < N^ is essentially satisfied in these systems. 



Bose-Einstein condensation of interacting atoms 



300 



200 



100 




100 200 300 

N (thousands) 



Fig. 2. Lack of saturation of the thermal component in a quantum degenerate 
atomic Bose gas. A'^' (black points) and A^o (grey points) are plotted versus the total 
atom number A*' at T = 177 nK and a = 135 oq. The corresponding predictions for 
a saturated gas are shown by black and grey solid lines. The critical point N — Nc 
is marked by a vertical dashed line. (Figure adapted from [6].) 



with just the value of the bound on the right-hand side slightly modified by 
interactions. However, as shown in Fig. 2, this is far from being the case under 
typical conditions of an ultracold gas experiment. Here, in an experimental 
series taken with a — 135 ao (where ao is the Bohr radius) and T — 177 nK, 
the measured critical atom number is N^ ~ 200 000; if the total number of 
atoms is increased to 450 000, only half of the additional atoms accumulate in 
the condensate. 

In order to explore the relationship between the experimentally observed 
non-saturation of the thermal component and the interatomic interactions, 
we first identify the relevant interaction energy. As a BEC is formed and then 
grows, the change in the average density of the condensed atoms is much 
larger than the change of the thermal density. Therefore one expects the non- 
saturation of the thermal component to result primarily from its interaction 
with the condensate (the 2gna{r) term in Eq. (10)). The relevant energy scale 
is then provided by [9] 



Mo ^ 9no{r = 0) = 



huj 



15iVo- 



fllK 



2/5 



(11) 



where aho = (h/mujy^^ is the spatial extension of the ground state of the 
harmonic oscillator. The energy /iq is the MF result for the chemical potential 
of a gas with iVo atoms at zero temperature in the Thomas- Fermi limit [9] . 
Guided by this scaling, for the data shown in Fig. 2 we plot N' as a 



.2/5 



function of A^'q' in Fig. 3. The growth of A^' with Nq' is not perfectly linear, 
so we quantify the non-saturation effect with two linear slopes: (1) the initial 
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Fig. 3. Quantifying the lack of saturation. Here A'^' is plotted as a function of 



N, 



g for the same series as in Fig. 



2. The horizontal dotted hne is the saturation 



prediction A''' = A^c- The two black lines show the initial slope So and the slope 
S for 0.1 < Ho / k^T < 0.3. The solid grey line is a guide to the eye based on a 
second-order polynomial fit. (Figure adapted from [6].) 



slope 5*0 for Nq -^ 0, and (2) the course grained slope S = A[N']/A[Nq''^] for 
0.1 < /io / I^bT < 0.3 [6]. The data shown in Fig. 3 can also be described by a 
second-order polynomial fit, as described later. 

The initial slope So may be compared with the HF model. In order to 
obtain the non-saturation effect to first order within the HF approach, we only 
consider the repulsive interaction of the thermal atoms with the condensate 
and not with other thermal atoms. From Eq. (10) this leads to an effective 
potential''' 

Feff(r) = T/(r) -h 2gno{r) = \V{r) - ^io\ + A^o • (12) 

Note that within this theory Nc = N^ since 14ff(r) = V{r) when Nq ~ 0. By 
integrating Eq. (7) with the effective potential of Eq. (12) one can predict a 
linear variation of N' /N^ with the small parameter /io / k^T: 



N' 
NO 



= 1 



Mo 
knT 



(13) 



with a ~ C(2)/C(3) « 1.37. This first order non-saturation result is identical 
to that obtained in more elaborate MF approximations, which only modify 
higher order terms. 

The origin of the non-saturation effect can be qualitatively understood by 
noting that interactions with the condensate modify the effective potential 
seen by the thermal atoms from a parabola into the "Mexican hat" shape 
of Eq. (12); this effectively allows the thermal component to occupy a larger 
volume, which grows with increasing A^q- 



Note that gno(r) = max{/io — V{r), 0}. 
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Fig. 4. Deviation from the saturation picture at a range of interaction strengths 
and temperatures. The non-saturation slope S is plotted versus the dimensionless 
interaction parameter X ex T'^a?'^ (see text). A linear fit (black line) gives dS/dX = 
2.6 ± 0.3 and an intercept 5(0) = —20 ± 100, consistent with complete saturation 
in the ideal-gas limit. The data points are based on measurements with the "^^K gas 
(closed circles) at a range of scattering lengths {a — 40 — 356 ao) and temperatures 
(T = 115 — 284 nK), and two additional experimental series taken with a *^Rb gas 
(open circles). (Figure adapted from [6].) 



From Eqs. (11) and (13) wc define the HF non-saturation slope 

dN' C(2) 



Si 



HF 



diK 



2/5 X 



C(3) 



X 



where X is the dimensionless parameter 
X 



C(3) fkBTV /15a^^/^ 



hbJ 



a\u 



(14) 



(15) 



The measured ^o is found to agree with 5'hf for a range of a and T values [6] . 

We now consider the non-saturation at higher A^o values, where the data 
is not well described by Eq. (14). Fig. 4 summarises the non-saturation slopes 
S{a, T) for a wide range of interaction strengths and temperatures. Within 
experimental error all data points fall onto a straight line with gradient 2.6 ± 
0.3 and intercept S{0) = —20 ± 100 when plotted against the dimensionless 
interaction parameter X. 

The first and most important thing to notice is that both non-saturation 
slopes, 5*0 and S, tend to zero for X — > 0. These experiments thus confirm the 
concept of a saturated Bose gas, and Bose-Einstein condensation as a purely 
statistical phase transition in the non- interacting limit. 

A question that is still open is the deviation of S from the first order 
HF resuh S'hf, i.e. dS/dX « 2.6 versus C(2)/C(3) « 1-37 predicted from HF 
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theory. This discrepancy can partially be explained by higher order terms in 
the mean- field theory, either directly from using Eq. (12) in Eq. (7), or using 
more elaborate MP theories such as the Popov approximation [3] . However the 
effect is far stronger in the experimental results than any of these MP theories 
predict. To see this we consider the next order term in Eq. (13), writing 
N'/N^} = 1 + ai^io/knT) + a2ifXo/kBT)'^ with a = 1.37. Experimentally, 
by identifying 5* with the gradient of this quadratic function evaluated at 
IJ-o/kBT = 0.2, we get a2 = 3 ±0.7. For comparison the Popov approximation 
gives a2 ~ 0.6. At present the reason for this discrepancy and the possible 
role of beyond-MP effects are unclear, and require further investigation. 

In summary, we can quantify the non-saturation of the thermal component 
in a harmonically trapped gas by writing the number of thermal atoms in a 
partially condensed cloud as 

N'^N, + SoN^/^ + S2N^'^ , (16) 

where 



_ C(2) fkBTy fl5a 
and 



2 \ hu J \ahoJ 

We have seen that in a harmonic trap the dominant non-saturation effect 
is "geometric" , arising from an interplay of the mean-field repulsion and the 
inhomogeneity of the condensate density. It is then interesting to consider the 
case of a uniform system, where this geometric effect is absent. Within MP 
theory, as the total density of a uniform Bose gas is increased past the critical 
value the thermal density n' actually decreases. This is due to the fact that the 
atoms in the condensate have less interaction energy, as discussed in section 
1.2. In addition, close to the transition, beyond-MP effects are expected to 
play an important role. This would make measurements in a uniform system 
(or of the local density in a trapped system) particularly interesting. 



4 Interaction shift of the transition temperature 

Having considered the effect of interactions on the saturation of the thermal 
component we now consider the location of the critical point itself. 

It is generally accepted that in a uniform system there is no interaction 
shift of the critical temperature Tc at the level of mean-field theory. However, 
the beyond-MP correlations between particles which develop near the critical 
point are expected to shift T^ [10,17-28]. Por several decades there was no 
consensus on the functional form, or even on the sign of this Tc shift (for an 
overview see e.g. [10,23,26,28]). It is now generally believed that the shift is 
positive and to leading order given by [23,24]: 
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Fig. 5. Opposing effects of interactions on tlie critical point of a trapped Bose gas. 
We sketch the density distribution in a harmonic potential V{r) at the condensation 
point. Compared to an ideal gas (dotted line) at the same temperature, repulsive 
interactions reduce the critical density, but also broaden the density distribution 
(solid line). Mean-field theory (dashed line) captures only the latter effect, and 
predicts an increase of the critical atom number Nc at fixed T, equivalent to a 
decrease of Tc at fixed N. (Figure adapted from [7].) 
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where AT, — T, — T^ and Ao is the thermal wavelength at temperature T^. 
Equivalently, the n^ shift at constant T is Auc/n"^ « -(3/2)Z\Tc/T°. The 
positive ATc implies that condensation occurs at a phase space density below 
the ideal-gas critical value of 2.612. 

The problem of the T, shift in a harmonically trapped gas is even more 
complex. In this case, at least for weak interactions, the shift is dominated by 
an opposing effect that reduces the critical temperature [29] . This negative T, 
shift is due to the broadening of the density distribution by repulsive interac- 
tions (see Fig. 5). To leading order it can be calculated analytically using MF 
theory [29], by self-consistcntly solving Eq. (6) with Ves = V^(r) + 2gn'{r): 



AT, 

TO 



-3.426 



An 



(20) 



For the experimentally relevant range of interaction strengths, < a/Ag < 
0.05, we numerically obtain the second-order MF shift, w 11.7 (a/Ao)^ [8]. 

The two opposing effects of repulsive interactions on the critical point of 
a trapped gas arc visually summarised in Fig. 5, where we sketch the density 
distribution at the condensation point for an ideal and an interacting gas 
at the same temperature^. In the spirit of LDA, the critical density should 
be reduced by interactions. However, interactions also broaden the density 
distribution. For weak interactions the latter effect is dominant, making the 
overall interaction shift ANc{T) positive, or equivalently ATc{N) negative. 

The dominance of the negative AiF shift of T, over the positive beyond-MF 
one goes beyond the difference in the numerical factors in Eqs. (19) and (20). 



^ The shift of the critical point can be equivalently expressed as ATc{N) or ANc{T), 
with ANciT)/N° « ~3An/T°. 
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In a harmonic trap, at the condensation point only the central region of the 
cloud is close to criticality^; this reduces the net effect of critical correlations 
so that they affect Tc only at a higher order in a/Ao. The MF result of Eq. (20) 
should therefore be exact at first order in o/Aq. The higher-order beyond-MF 
shift is still expected to be positive, but the theoretical consensus on its value 
has not been reached [30-34]. 

4.1 Measurements on a harmonically trapped Bose gas 

Since the early days of atomic BECs there have been several measurements 
of the interaction T^ shift in a harmonically trapped gas [35-37]. These ex- 
periments, performed at a/Ao ranging from 0.007 [37] to 0.024 [36], nicely 
confirmed the theoretical prediction for the linear MF shift of Eq. (20), but 
could not discern the beyond-MF effects of critical correlations. 

The recent measurements [7] presented here provided the first clear ob- 
servation of the beyond-MF Tc shift in a trapped atomic gas. Several im- 
provements contributed to making this possible. First, we explored slightly 
higher interaction strengths, up to a/Ao ~ 0.04. Second, by performing pre- 
cision measurements outlined in section 2, and directly accessing the small 
differential Tc shift due to the variation in a/Ag, we significantly reduced the 
experimental error bars. Third, understanding the non-saturation effects dis- 
cussed in section 3 was also essential for accurately determining the critical 
point from the measurements performed close to it (see Fig. 6). 

Fig. 6 illustrates a differential measurement with a = 56 oq and a — 27 A a^. 
The rise of Nq in Fig. 6(a) is not simply vertical because the thermal compo- 
nent is not saturated at Nc (sec section 3). It is therefore essential to carefully 
extrapolate N' to the A^o = limit in order to accurately determine Nc- 
For small iVg the extrapolation is done using N' = Nc + SqNJ , with the 
non-saturation slope S'o(T', w,a) given by Eq. (17). 

In Fig. 7 we summarise our measurements of ATc/Tc [7]. The data taken 
with different atom numbers, N k, (2 — 8) x 10^, fall onto the same curve, 
confirming that the results depend only on the interaction parameter a/Ao. 
The MF prediction agrees very well with the data for a/Ao ^ 0.01, but for 
larger a/Ao there is a clear deviation from this prediction. All the data are 
fitted well by a second-order polynomial 

with 6i = —3.5 ± 0.3 and 62 = 46 ± 5. The value of 61 is in agreement with 
the MF prediction of —3.426. The &2 value strongly excludes the MF result of 
6^^ PS 11.7 and its sign is consistent with the expected effect of beyond-MF 
critical correlations. Fixing 61 = —3.426 (which is expected to be exact even 
including beyond-MF effects) gives an improved estimate of 62 = 42 ± 2. 



^ The size of the central critical region is Tc ~ (a/Ao)i?T, where Rt ~ ^/ksT/r 
is the thermal radius of the cloud [30] . 
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Fig. 6. Determination of the critical point and the differential interaction shift, (a) 
Condensed (A'^o) versus thermal (A'^') atom number for two concurrently taken data 
series with a — 56aQ (circles) and a — 274 ao (squares). Note that all points corre- 
spond to condensed fractions below 2 %. The data is scaled to the same temperature 
(T = 240 nK) and shows the shift of the critical point in the form ANc{T). Solid 
lines show the extrapolation to A'o = 0, necessary to accurately determine Nc- (b) 
A'^' is plotted versus A^q for the same data as in (a), showing more clearly the 
extrapolation procedure. (Figure adapted from [7].) 
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Fig. 7. Interaction shift of the critical temperature. Data points were taken with 
N ~ 2x IC (open circles), 4x 10"" (black squares), and 8x 10' (open triangles) atoms. 
The dashed line is the MF prediction. The solid line shows a second-order polynomial 
fit to the data (see text). Vertical error bars are statistical, while systematic errors 
in ATc/T^ are assessed to be < 1 % [7]. (Figure adapted from [7].) 



4.2 Connection with a uniform Bose gas 

In order to make a connection between the experiments on trapped atomic 
clouds and the theory of a uniform Bose gas we also need to consider the effect 
of interactions on the critical chemical potential fXc- 

In a uniform gas the interactions differently affect Tc (or equivalently Uc) 
and /ic at both MF and beyond-MF level. The simple MF shift /3^^^ = 
4 C(3/2) a/Ao (where /3 = l/k^T) has no effect on condensation. To lowest 
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n(x) 



~ a/lo 




Fig. 8. Beyond-MF effects near tlie critical point in a harmonically trapped Bose 
gas. (a) For a fixed T, the density distribution at the critical point N = Nc < N^^ 
(solid grey line) is compared with the MF prediction (dashed line). In the trap 
centre we expect n^''^ — ric oc a/Ao, characteristic of a uniform system. However the 
experimentally measured A*'^^ — Nc cc (a/Ao)^ is dominated by the density shift 
outside the central critical region, and is not directly related to the ric shift, (b) If 
N is increased to A^^'' > Nc, a condensate induced by critical correlations forms 
within the critical region of size oc a/Xo- The condensed atom number A^o oc (a/Ao)'' 
directly relates to the critical density shift Auc oc a/Ag. (Figure adapted from [8].) 



beyond-MF order we have^: 



I3fi, « /^A^r* + B2 



Ar 



(22) 



We see that there is a qualitative difference between Eqs. (19) and (22). Specif- 
ically, we have n^^ — ric (x a/Xo, but /x^^ — /Xc oc (a/Ao)^. This difference 
highlights the fact that the problem of the Tc shift is non-perturbative; near 
criticality the equation of state does not have a regular expansion^ in /x, oth- 
erwise one would get Auc ex /ic — M^^- 

For a harmonic trap, within LDA the uniform-system results for n^ and 
He apply in the centre of the trap, and elsewhere local /i is given by Eq. (5). 
The result for the Tc shift however does not carry over easily to the non- 
uniform case. As illustrated in Fig. 8(a), in the centre of the trap we expect 
Arte oc a/Xo, but the measured beyond-MF ANc oc (a/Ao)^ (see Eq. (21)). 

In fact, the experimentally observed T^ shift (Eq. (21)) qualitatively mir- 
rors Eq. (22), with the MF term linear in a/Ag and the beyond-MF one 
quadratic in a/Ao (to leading orders). This similarity can be understood as 
follows: The interaction shift of fic affects the density everywhere in the trap; 
outside the small critical region the equation of state is regular in /i and the 
local density shift is simply proportional to the local fi shift; the contribution 

* Note that B2 is not just a constant but contains logarithmic corrections in a/Ao 
[23]. We neglect these in our discussion since they are not discernible at the 
current level of experimental precision. 

^ The non-interacting equation of state (Eq. (3)) cannot be expanded about Dc 
in Pfi, but rather in ^—jSii; up to first order this expansion gives D — Dc — 
2\/^V^/^M- This scaling goes some way in explaining the qualitative difference 
between Eqs. (19) and (22) although it cannot be used quantitatively. 
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to the total Nc shift from the non-critical region greatly outweighs the contri- 
bution from within the critical region and therefore one qualitatively expects 
Nc — N^^ (X ^c~ M^^- More quantitatively, this connection is given by^° [23] 

-3(&2-6r) = ^S2. (23) 

The key conclusion of this discussion is therefore that the bcyond-MF Tc 
shift observed in a trapped gas is directly related to the beyond-MF fic shift 
(in either trapped or uniform system). It does not however directly reveal the 
expected linear n^ shift and the theoretically most intriguing non-perturbative 
connection between fic and Uc shifts. In fact one could say that the historically 
puzzling (theoretical) connection between ric and iic is just replaced by the 
puzzling connection between (expected) Uc and (measured) Nc- 

This problem can be partly overcome by studying the condensed fraction 
/o in a trapped atomic cloud at the MF-predicted critical point [8]. By def- 
inition /o vanishes within MF theory, and so directly measures the effect of 
critical correlations which shift Tc above T^^ . At a fixed T we consider the 
condensed fraction of a gas, Nq/N, at the point where N = iV^'^ > Nc, 
as illustrated in Fig. 8(b). The analogous quantity for a uniform gas was 
first theoretically studied by Holzmann and Baym [38], who showed that^^ 
na/n ex Ahc oc o/Aq. The condensate density tiq vanishes at the point where 
the local n = fic, so from Eqs. (5) and (22) we get that the spatial extension 
of the condensate is also ex a/Ao, and hence /o oc (a/Ao)"*. 

The quartic scaling of /o with a/Ao is thus directly related to the expected 
linear scaling of Auc with a/Ag; the two scalings are simply connected by 
the volume of the critical region, oc (a/Ao)"^. This scaling is indeed confirmed 
experimentally [8], as shown in Fig. 9. 

One can take the comparison with theory beyond just the scaling of /o 
with a/Ao, and quantitatively compare the measured A'o at the MF critical 
point with the Monte-Carlo (MC) calculations [39] for a uniform gas. This also 
works out very well, with the measured and predicted Nq agreeing within a 
few percent [8] . It is however important to carefully summarise the conceptual 
steps involved in this comparison: 

(i) On the one hand, the measurements of the Tc shift [7] experimentally 
provide (up to logarithmic corrections) the value of /x^^ — /Zc oc (a/Ao)^, via 
Eqs. (21), (22), and (23) [8]. 

(ii) On the other hand, the MC calculations [24,39] that predict the Uc shift 
of Eq. (19) also provide tabulated values of the uniform-system condensate 
density hq for any /i — /ic oc (a/Ag)^. 



^° Note that this relationship is closely related to Eq. (13). 

^^ This scaling holds for any distance from the critical point given by (p — 
Mc)(Ao/a)^ ~ const. By applying it to the MF critical point we neglect the loga- 
rithmic corrections to /i^^ — fic, which are so far not experimentally observable. 
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Fig. 9. Condensed fraction of an atomic gas induced by critical correlations. The 
condensed fraction /o = Nq/N is measured for A'^ = A^^'' > Nc- A fit to the data 
(solid line) with the function /o oc (a/Ao)^ gives an exponent x — 3.9 ± 0.4, in 
agreement with the predicted x — 4. (Figure adapted from [8].) 



(ill) Combining these two results and the LDA (Eq. (5)) we calculate the 
expected Nq in a trapped gas with N ~ N^^ and find excellent agreement 
with the measurements shown in Fig. 9 [8]. 

Overall this provides strong cross-validation of theory and experiment. 
However it is important to note that the two different measurements of 
beyond-MF effects (i.e. of the T^ shift and /o) do not provide two independent 
quantitative tests of the uniform-system theory. Instead, what we have shown 
is that they are consistently connected via the MC calculations for a uniform 
system. Finally it is also important to stress that we are still lacking a direct 
measurement of the ric shift, which would explicitly test the historically most 
debated theoretical result of Eq. (19). This goal remains open for future mea- 
surements, either of the local density in a harmonically trapped gas or on a 
uniformly trapped atomic gas. 



5 Equilibrium criteria and non-equilibrium effects 

Finally, we discuss the criteria for the measurements on a trapped Bose gas to 
faithfully represent its equilibrium properties, and the non-equilibrium effects 
revealed when they are violated. It is helpful to distinguish two types of non- 
equilibrium behaviour, transient and intrinsic. 

Transient non-equilibrium effects are more familiar, and occur whenever 
some system parameter, such as the interaction strength, is rapidly changed 
(quenched). After such a quench the system relaxes towards its new equilib- 
rium. Classically, one can estimate the relaxation time to be several elastic 
scattering times l/7ei, where 7oi is the elastic scattering rate [40-43]. 

However a system with continuous dissipation can only be "close to" ther- 
modynamic equilibrium and is to some extent always intrinsically out of cqui- 
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Fig. 10. Non-equilibrium effects, (a) ATc/T^ is determined assuming equilibrium, 
as in Fig. 6. At both very low and very high a the apparent Tc deviates from 
the equilibrium curve, (b) Equilibrium criteria (see text); 7cir (solid squares) is 
the number of elastic collisions per particle during 1 % atom-loss; -^Ici/Oj = 1 (open 
circles) marks the onset of the hydrodynamic regime. (Figure adapted from [7].) 



librium. The proximity to equilibrium broadly depends on the competition 
between relaxation and dissipation. For an atomic gas, this leads to a crite- 
rion based on the dimensionlcss parameter 7oit, where t is some characteristic 
dissipation time, e.g. for atom loss. In practice, the relevant t and the criteria 
for equilibrium measurements depend on the required measurement precision. 

In the case of the Tc measurements presented in section 4, Nc is determined 
to about 1 %, so we require that the gas continuously (re-)equilibrates on a 
timescale t corresponding to only 1 % atom-loss. This requires about 100 times 
higher 7ci than one would naively conclude by taking the 1/e lifetime of the 
cloud as the relevant timescale. 

An interesting question is then what happens if we violate these stringent 
equilibrium criteria. In Fig. 10(a) we show measurements extending beyond 
the equilibrium region shown in Fig. 7, and in Fig. 10(b) we plot the cor- 
responding 7ciT. Individually, 7ci and t vary vastly as a function of a (7ci 
increasing and r decreasing) [7] , but the breakdown of equilibrium appears to 
occur at 7cit « 5 in both the low- and high-a limit. 

In the small-a limit the apparent Tc is significantly above the equilibrium 
curve. We can qualitatively understand this effect within a simple picture. In 
this regime, losses are dominated by one-body processes which equally affect 
A'o and N' . The net effect of equilibrating elastic collisions would therefore 
be to transfer atoms from the condensate to the thermal cloud. However the 
dissipation rate is too high compared to 7oi, and so A^o remains non-zero even 
after the total atom number drops below the equilibrium critical value Nc- 



18 Robert P. Smith and Zoran Hadzibabic 

In the large-a limit the initial breakdown of equilibrium again appears to 
result in condensates surviving above the equilibrium T^. However the physics 
in this regime is much richer, with several potentially competing effects re- 
quiring further investigation. For example, three-body decay affects A^o and 
A'^' differently, the thermal component is far from saturation, and the thermal 
component of the gas also enters the hydrodynamic regime, ^c\/uj > 1. 

It is interesting that we observe "superheated" BECs for both very weak 
and very strong interactions. Further study of these effects should prove use- 
ful for improving the understanding of condensation in intrinsically out-of- 
equilibrium systems, such as polariton gases. 
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